Proximity effect, quasiparticle transport, and local magnetic moment in 
ferromagnet-d-wave superconductor junctions 
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The proximity effect, quasiparticle transport, and local magnetic moment in ferromagnet-d-wave 
superconductor junctions with {110}-oriented interface are studied by solving self-consistently the 
Bogoliubov-de Gennes equations within an extended Hubbard model. It is found that the proximity 
induced order parameter oscillates in the ferromagnetic region. The modulation period is shortened 
with the increased exchange field while the oscillation amplitude is depressed by the interfacial 
scattering. With the determined superconducting energy gap, a transfer matrix method is proposed 
to compute the subgap conductance within a scattering approach. Many novel features including the 
zero-bias conductance dip and splitting are exhibited with appropriate values of the exchange field 
and interfacial scattering strength. The conductance spectrum can be infiuenced seriously by the 
spin-fiip interfacial scattering. In addition, a sizable local magnetic moment near the {110}-oriented 
surface of the d-wave superconductor is discussed. 



PACS numbers: 74.20.Mn, 74.80.Fp, 74.50. +r 
I. INTRODUCTION 

The electronic trans- 

port in ferromagnetic-superconducting hybrid structures 
is currently a very active area of research due to their 
interesting physical properties and potential device ap- 
plications. A fundamental transport process through the 
interface between the normal conducting and supe|E|Con- 
ducting materials is the Andreev reflection (AR):EJ An 
electron incident with energy below the superconduct- 
ing energy gap cannot enter the superconductor, it is in- 
stead reflected at the interface as a hole by transferring a 
Cooper pair into the superconductor. The earlier spin po- 
larization experiments involving ferromagnet and. super- 
conductor were performed on tunnel junctionsp where 
the AR is unimportant due to the strong interface bar- 
rier. Recently, the effect of spin polarization on the AR 
has been investigated in ferromagnet^jCpnvcntional su- 
perconductor contacts experimentally,BQ where the AE. 
plays an important role. In an earlier theoretical work,El 
this effect was studied in the zejjOsjbias limit. Several re- 
cent spin injection experimentsa El have been done with 
high-Tc superconductors. Common to both ferromagnet- 
conventional and high- Tc superconductor junctions, the 
subgap conductance at a given bias is suppressed due 
to the suppression of AR by the spin splitting of energy 
band in the ferromagnet. In particular, a zero-bias con- 
ductance dip was obserV|ed in the ferromagnet-high-Tc 
superconductor contacts.la It has been widely accepted 
that the high-Tc superconductors have a d-wave pairing 
symmetrylij The above interesting observation may indi- 
cate the importance to take into account the unconven- 
tional pairing symmetry of the cuprate superconductors. 
It has been shown that, due to the fcppation of midgap 
states at the {110}-oriented interface Jld the conductance 
spectrum of d-wave superconductor junctions differs dra- 
matically fro ni t ha t of conventional s-wave superconduc- 
tor junctions.lla't2l Thus the difference should also exist 



between ferromagnet-d-wave and s-waiffijSuperconductor 
junctions. In recent theoretical works,li3t£l the novel fea- 
tures of AR have been exhibited in the subgap conduc- 
tance of ferromagnet-d-wave superconductor junctions. 
More recently, the effect of spin injection into s- and d- 
wave superconductors has also been studied with an em- 
phasis on the interplay between boundary and bulk spin 
transport processes. Ell 

In parallel, there also has been much interest in the 
interplay of superconductivity and ferromagnetism in 
these combined structures. In the case of ferromagnet- 
superconductor multilayers, the transition temperature 
of the (s-wave) superconductor changes nonmonotoni- 
cally with the thickness or the jexchange field strength 
of the ferromagnetic layers .EaEil In superconductor- 
ferromagnet-superconductor junctions, the exchange 
field in the ferromagnetic layer lieads to oscillations 
of the Josephson critical current El EJ More recently, 
the influence of the exchange field on the Josephson 
current in superconductor-ferromagnet-superconductor 
junctions with uncttnventional pairing symmetry has 
also been studied.Ej In ferromagnet-superconductor- 
ferromagnet multilayers, the appearance of the super- 
conducting energy gap causes a reduction of the indicept 
magnetic coupling which exists in the normal stateEj'EZI 
For the case of ferromagnet-superconductor junctions, 
the superconducting proximity effect will also change 
in thfer^rpaence of an exchange held. The previous 
workgallinla with an emphasis on the transport through 
the ferromagnetic-superconductor junctions were based 
on a simplified continuum model and did not calcu- 
late the order parameter (i.e., pairing amplitude) self- 
consistently so that the proximity effect cannot be in- 
cluded. 

The purpose of this work is to present a unified and 
rigorous treatment of the proximity effect, transport 
and magnetic properties in a ferromagnet-d-wave super- 
conductor junctions. Within the framework of an ex- 



tended Hubbard model, we solve self-consistently the 
Bogoliubov-de Gennes (BdG) equations to obtain the 
spatial variation of the order parameter and supercon- 
ducting energy gap. With the obtained energy gap, a 
transfer matrix method is then proposed to calculate the 
subgap conductance within the scattering approach. The 
self-consistent calculation also allows us to study the local 
magnetic moment in the superconducting region due to 
the presence of the exchange field or Zeeman coupling. 
The main procedure in the present work is parallel to 
the study of transport properties in the normal-metal- 
anisotropic superconductor junctions using the quasiclas- 
sical theory, where the pair potential first obtained from 
the quasiclassical formalism is substituted into the An- 
dreev equation to, calculate the reflection and transmis- 
sion coefficientsra 

The paper is organized as follows: In Sec. |l| the BdG 
equations for the ferromagnet-superconductor junctions 
are derived within the extended Hubbard model. In 



W = -t ^ cLcjo- + ^ UiTiia + ^ hi^Tiic, -/iy^i 



Sec. [H, the order parameter and pair potential are de- 



termined self-consistently. The subgap differential con- 
ductance and the local magnetic moment are presented 
in Sees. |^ and 0, respectively. Finally, conclusions are 
given in Sec. p|. 



II. THE BOGOLIUBOV-DE GENNES 

EQUATIONS FOR THE 

FERROMAGNET-SUPERCONDUCTOR 

JUNCTIONS 

We use a single-band extended Hubbard model to de- 
scribe the ferromagnet-superconductor junctions. The 
geometry is shown in Fig. m for (a) a ferromagnet-s- 
wave superconductor junction with a {100}-oriented in- 
terface and (b) a ferromagnet-(i2,2_ 3,2 -wave superconduc- 
tor junction with a {110}-oriented interface. The signifi- 
cant difference between s-wave and d-wave superconduc- 
tors can be exhibited most clearly in these two struc- 
tures. For such a crystalline orientation of the (1^2 _^2- 
wave superconductor, a dxy-wa.ve superconductor junc- 
tion is formed. For the junction involving the s-wave 
superconductor, the qualitative features are indepen- 
dent of the interface orientation. Hereafter, we will call 
the ferromagnet-s-wave superconductor junction with 
the {100}-oriented interface the FS junction while the 
fcrromagnet-(i2.2_2.2-wave superconductor junction with 
the {110}-oriented interface the FDa;j, junction. In the 
junction geometry, both the ferromagnet and the super- 
conductor are treated as semi-infinite. Here we choose 
the interface to be at the 0-th layer. We further assume 
that the transition temperature of the superconductor is 
much smaller than the Curie temperature of the ferro- 
magnet so that fluctuation effects on the magnetism are 
negligible. 

Under these assumptions, the Hamiltonian defined on 
two-dimensional square lattice is given by 



(ij)a- icr icT ia 

-^M^)ni]nii - - Y^ Vi{ij)ni^nj„> . (2.1) 

i (ij)o-CT' 

Here i and j are site indices and the angle bracket im- 
plies that the hopping and interactions are only consid- 
ered up to nearest-neighbor sites, cJ^ {cia) are creation 
(annihilation) operators of an electron with spin a on 

site i, nio- = Cj^Cio- is the electron number operator on 
site i, t the hopping integral, and fi the chemical poten- 
tial. The interfacial scattering potential is modeled by 
Ui — UoSno^ where n is the layer index along the direc- 
tion perpendicular to the interface plane. The conduction 
electrons in the ferromagnet interact with an exchange 
field, hicr = —hQazQ{—n), where 8 is the Heaviside step 
function and cr^ {— ±1) is the eigenvalue of the z com- 
ponent of the Pauli matrix. A real space representation 
of the exchange interaction is used since the present sys- 
tem is inhomogeneous. The quantities Vb(i) and Vi(ij) 
are on-site and nearest-neighbor interaction strength, re- 
spectively. They are taken to be Vq and Vi in the super- 
conductor and identically zero in the ferromagnet. Pos- 
itive values of Vq and Vi mean attractive interactions 
and negative values mean repulsive interactions. When 
Vb < and Vi > 0, the d-wave pairing state is favorable. 
Here we also would like to point out that by taking the 
same chemical potential in both the ferromagnet and the 
superconductor, we have ignored the effect of the Fermi 
wavevector mismatch between two materials. Very re- 
cently, this effect on the conductance spectrum has been 
well studied within the simple continuum model.E£l 

Within the m ean- field approximation, the effective 
Hamiltonian Eq. (|2.l|) can be written as 



Weft = -i'Y 4o-'^Jo- + X! ^''^■'^ + X! ^i'^'^i'^ -^y^nia- 

(ij)cr io- i(T i(T 

+ I][AS(i)cij,Cit-f Ao(i)ct^ctj 

i 

+ Y. [Ai(i)cjiCiT-KA^(i)c|^c]j, (2.2) 



where 



Ao(i) = K)(i)(ciTCii 



A^(i) = T/i(i,i + J)(ciTCi^^ 



+0i/ 



(2.3) 



(2.4) 



are the on-site and nearest-neighbor pair p oten tials, re- 
spectively. The effective Hamiltonian Eq. (2.2) can be 
diagonalized as H^k = Eg + J2u Eiultlu by performing 
the Bogoliubov transformation. 
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(2.5) 



Here 7J (7;^) are creation (annihilation) operators of 
fermionic quasiparticles. E^, are the quasiparticle eigen- 
values. The quasiparticle wavefunction-, amplitudes 
{uiiT,Via) satisfy the lattice BdG equationsEj 
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(2.6a) 



(2.6b) 



(2.7) 



direction. The problem becomes solving the BdG equa- 
tions for {un{ky), Vn{ky)) corresponding to the eigenvalue 
E{ky): 
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(2.12a) 



(2.12b) 



where 

Hnn' = -2t COs{kya/ V2)5n±l^n' + {Un ~ fJ-)Snn' , (2.13) 



Aij = Ao(i)(5ij+A^(i)5.^j_ 



(2.8) 



with S = ±Xa, ibxft the unit vectors along the crystalline 
Xa and Xb axis. The energy gaps for on-site and nearest- 
neighbor pairing are determined self-consistently 



Ao(i) 



Vo(i)i^o(i) 



[<.<!* + u^iT] tanh(^,/2T) , (2.9) 



■ /^["-iT^n 



Aj(i) 



^i(i,i- 
V^i(i,i- 



J)F^(i) 



2 Z^^"iT^i+J,i 

X tanh(£'^/2r) , 



^+5a^iT 



(2.10) 



where the Boltzmann constant ks — I has been taken, 
and Fq (i) and Fg (i) arc the on-site and nearest-neighbor 
bond order parameter. Note that the 4x4 BdG equations 
are decoupled into two sets of 2 x 2 equ atio ns since the 
spin- flip effect is not considered in Eq. (2J). Note that 
the eigenstates of the BdG eq uatio ns exist in pairs: // 
{ui,vi) is the solution of Eq. ( 2.6a ) with the eigenvalue 
E, then {—v't,ut) is the solution of Eq. (276a) with the 
eigenvalue —E. 

For a clean ferromagnet-superconductor junction with 
a flat interface, which we are considering, there exists the 
translation symmetry along the y direction so that the 
Bloch theorem can be applied to this direction. Then for 
the FDxy junction, the eigenfunction can be written in 
the form 






/K, 



Un(Ky) 



Jimkya I \/2 



(2.11) 



where a is the lattice constant, ky e \—tt j \f2a^'K j \f2a\^ 
m and n are ionic-layer indices in the x and y direc- 
tions, and iVj, is the number of unit cells along the y 



A„„, = Ao(n)5„„, + [Aa(n,n± 1)6^''="''/^ 

+Afc(n, n ± l)e±^^-''''/^]5„±i,„. , (2.14) 

with the gap functions given by 

Vo{n) 



Ao(n) 



2iV„ 



^Kt<I + <i<T] tanh[£;,(fc,)/2T] 



u.ky 



Aa(n,n±l) 



(2.15) 



Vi{n,n±l) ^ ±.fc„a/v^ 

2N Z^KT^n±i,ie 



„a/V2l 

xtanh[£:^(fcj,)/2r] , 






(2.16) 



and 

Afa(n,n±l) 



Vi{n,n±l) y-^ ^,kya/^ 
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xtanh[£;^(fcy)/2T] . 



(2.17) 



The problem with other orientations of the flat interface 
can be treated similarly. 



III. SELF DETERMINATION OF THE ORDER 
PARAMETER AND THE PAIR POTENTIALS 

We solve the BdG equations self-consistently by start- 
ing with an initial gap function. After exactly diago- 
nalizing Eq. (2.6), the obtained Bogoliubov amplitudes 
are substituted into Eqs. (2.9) and (2.10) to compute a 
new gap function. We then use it as input to repeat the 



above process until the relative error in the gap func- 
tion between successive iterations is less than the desired 
accuracy. Throughout this work, we concentrate on the 
zero temperature case unless specified explicitly, and take 
the parameters: n = 0, and Vi — 2i and Vq — ~2t 
for d-wave superconductors while V^i = and Vq = 2t 
for s-wave superconductors. This set of parameter val- 
ues give the zero-temperature energy gap A^o = 0.241i 
and Aso = 0.376i for the bulk d-wave and s-wave super- 
conductors, respectively. The corresponding cohereaee 
$,d » 1.3a and ^s « 3.4a. Note that as in other works,Ej'EJ 
the model parameters chosen here are not intended for 
realistic materials. 

For the d-wave superconductor, the amplitudes of d- 
and extended s-wave order parametecs can be defined in 
terms of the bond order parameters:E3 

Pd{i) = J[^s.(i) + i^-x„(i) - F,,{i) - F_£,(i)] , (3.1a) 
^.(i) - ^[^^..(i) + F-iA^} + F,,{i) + F^,,{i)] . (3.1b) 

Accordingly, the d-wave and extended s-wave pair poten- 
tials are given by: 

Ad(i) = i[A£„(i) + A_s„(i) - As,(i) - A_i,(i)] , (3.2a) 

A,(i) - i[As„(i) + A_,„(i) + As,(i) + A_,,(i)] . (3.2b) 

In the superconducting region, the energy gap is propor- 
tional to the order parameter because of the constant 
pairing interaction. In the bulk state of d-wave super- 
conductor, the extended s-wave component is zero. For 
the junction systems under consideration, the induced 
extended s-wave component near the interface is numeri- 
cally found to be vanishingly small for the value of on-site 
repulsive interaction we have taken. For the conventional 
s-wave superconductor, the order parameter and the pair 
potential are directly on-site defined. 

In Figures || and g, we plot the spatial variation of 
the order parameter for various values of exchange field 
in the FS junction and the FD^jj, junction. In this case, 
there is no interfacial scattering potential. As can be 
seen, the exchange field does not influence the order pa- 
rameter in the superconducting region. In the normal 
metal case {ho = 0), the proximity induced order pa- 
rameter monotonically decays into the normal metal re- 
gion. Interestingly, common to both the FS and FDxy 
junctions, the order parameter in the ferromagnetic re- 
gion no longer changes monotonically, it instead oscillates 
around the zero value of order parameter. In addition, 
as the exchange field is increased, the oscillation period 
becomes shorter. This behavior could be understood in 
the following way. Take the FS junction as an example. 
Since the component of the wavevector parallel to the 
interface is conserved, we can just consider the normal 
component. In the superconducting region, the wavevec- 
tors (momenta) of the spin-up and spin-down electrons 



forming the Cooper pairs have the same amplitude (but 
opposite directions) qx- Upon entering into the ferromag- 
netic region, the pair amplitude decays. Simultaneously, 
the spin- up electron lowers its energy by /iq, while the 
spin-down electron gain the energy Hq. In order for each 
electron to conserve its total energy, the spin-up electron 
should adjust its momentum to g^ , while the spin-down 
electron to qi- Therefore, from the expressions of the 
order parameter given by ( |2.9[ ), we can approximately 
write the order parameter as cos[(g| — qi)na]^{n), where 
$(«.) is a slow- varying envelope function. Therefore, the 
exchange field in the ferromagnet causes the spatial mod- 
ulation of the order parameter, which now roughly varies 
at the scale of {q^ — q^)^^. In the continuum model, the 
difference q] — qi ~ 2ho/hvFx, where vp^ is the nor- 
mal component of the Fermi velocity. This also explains 
the decrease of the modulation period with the exchange 
field ho. In Figs. and S, the spatial variation of order 
parameter are plotted for various values of the interfa- 
cial scattering potential in the FS and FD^y junctions 
with the exchange field fixed at ho = 0.125D {D = 8t 
is the band width). Our numerical results show that as 
the interfacial scattering potential becomes stronger, the 
oscillation amplitude of the order parameter in the fer- 
romagnet is decreased. This is because the amplitude of 
the slow-varying envelope function mentioned above has 
been suppressed by the interfacial scattering at the inter- 
face. However, in the superconducting region, the order 
parameters of the FS and FT)xy junctions show different 
behavior in the presence of the interfacial scattering. As 
shown in Figs. for the FS junction, the depression of 
the order parameter near the interface is decreased by the 
interfacial scattering. As the interface is strongly reflect- 
ing (large Uo), the superconductor and the ferromagnet 
are almost decoupled, and since the opaque interface it- 
self is not pair breaking for s-wave superconductivity, the 
s-wave order parameter is not depressed. In contrast to 
the s-wave case, in the FF)xy junction, the reflected quasi- 
particles from the {110}-oriented interface are subject to 
a sign change of the order parameter, which makes the 
opaque interface itself pair-breaking. Thus the d-wave or- 
der parameter is strongly depressed (see Fig. ||) . Finally, 
since we have assumed that there is no pairing interaction 
in the ferromagnet, the pair potential or energy gap in 
this region is zero. It is the pair potential that acts as an 
off-diagonal scattering potential in the BdG equations. 



IV. THE SUBGAP DIFFERENTIAL 
CONDUCTANCE 

A. The case of nonmagnetic interfacial scattering 



Once the BdG equations ( |2.6| ) are solved self- 
consistently, we can use the obtained pair potential to cal- 
culate the differential conductance. The transport prop- 
erties through the normal-metal-superconductor June- 



tions can be studied within the Blonder-Tinkham- 
Klapwijk (BTK) scattering formahsm,E3 which expresses 
the differential conductance in terms of the normal and 
Andreev reflection coefhcients. In contrast to the tun- 
neling Hamiltonian model, which requires an opaque 
barrier at the interface, the BTK theory can consider 
the case of an arbitrary barrier strength. Also no- 
ticeably, the BTK formalism can be regarded as the 
earliest version of the Landauer-BiittikeiH formula ap- 
plied to the coherent transport through a normal-metal- 
superconductor structure.Ej Recently, the BTK the- 
ory has been extended to the spin-dependent transport 
through ferromagnet-superconductor junctions^ Within 
the tight-binding model, the averaged differential con- 
ductance can be obtained as: 



G = 



hN„ ^^ 



[1 + Rh, 



Re 



(4.1) 



which shows clearly that an incoming electron of spin a 
(=1,1) is normally reflected as an electron of the same 
spin a with probability Re^^a = ko-crpj and Andreev re- 
flected as a hole of the opposite spin a with probability 
Rh.ajT = [siT^{qa_a)/ sin{qa-a)]\r^a-\'^- Here the summation 
is over all the transverse modes and over the spin in- 
dices. Without confusion, h in Eq. (4.1) is the Planck 
constant. In contrast to the continuum model, the fac- 
tor sm{qa_a) / sin{qcra) comes from the band structure ef- 
fect. Our previous work within the continuum modeltll 
concentrated on the direction-dependent subgap con- 
ductance through the ferromagnet-superconductor junc- 
tions, which can be experimentaily explored with the 
scanning-tunneling spectroscopy.Ej For a point contact 
junctionp a summation over the transverse modes is 
needed. 

The remaining thing is to obtain the Andreev and nor- 
mal reflection coefficients, which can be calculated using 
the transfer matrix method. As an illustration, we give a 
detailed procedure for the calculation of these coefficients 
for the FS junction, which has a {100}-oriented interface. 
From the BdG equations (p^) , we can write the relation 
of wavcfunctions among consecutive layers: 



(4.2) 



where the transfer matrix for n-th layer is given by 




M,„ 



/ 



V 



t t 

t t 

1 



1 



-1 \ 

-1 



0/ 



(4.3) 



with e„ = UoSno — 2tcoskya — fi and /i„ = hoQ{—n). 
The incident and transmitted wave amplitudes are then 
connected with the total transfer matrix 




Ml 




where 



MT = l[Mr 



(4.4) 



(4.5) 



Here Nl^u are the indices for two outmost layers of the 
scattering region where the energy gap has approached 
its bulk value. If the electron wave is incident with spin 
up and transverse momentum ky, the incident, reflected 
and transmitted wavcfunctions are written as 






UNl-I 
VNl-1 



■^iU 1 ) + nt 



(4.6a) 



J je-^'^^+ntf ? le-^'^'^ 



+nT 



o*9T° 



(4.6b) 



and 



--«^-tnf^+^^"M+^ /«-e^^ 



unr+i 

VNr + 1 



HT 



v+ 



in u^ ' ' ^^■'^''^ 



C|| 



v+ 



U^e^^ \^U.a^. IV^e 



HT 



^ — ikha 



(4.7b) 



Here i|| and i|| are the transmission amplitudes. (/)± are 
the internal phase of the energy gap and are identically 
zero for the s-wave superconductor. The wave vectors 
in the ferromagnet and the superconductor are, respec- 
tively. 



9Taa = cos" 
and 

ke.hO- = COS" 



2i cos kytt + ^± {E + ho) 
2t 



, (4.8) 



2t cos kya + fi± y/E'^ - |A(k±) 
2i 



(4.9) 



The BCS coherence factors are given by 



1 

2 


^ 1 V^2_|A(k±)p 
E 


(4.10a) 


1 
2 


V^2_|A(k±)|2 

E 


(4.10b) 



For clarity, we have written the energy gap explic- 
itly depending on the wavevector A(k±) = A{±qo,ky) 



with the a::-coinponent of the Fermi wavevector qoa = 
cos~^[—{2t cos kyU + n)/2t]. For a conventional s-wave 
superconductor, A(k±) = A^o- Then the reflection am- 
plitudes can be obtained by solving the linear equa- 
tion (4.4). The reflection coefficients for the FD^y 
junction can be calculated similarly, where qQa/y/2 ~ 
cos"i[-^/4tcos(fcj,a/\/2)] and A(k+) = -A(k_) = 

4Ado sin((j'oa/v2) sin(fcya/v2) so that the internal phase 
0+ = (or tt) while ^_ = tt (or 0) depending on A(k+) 
being positive or negative. Note that the reflection coeffi- 
cients, which we are calculating here, can also be used|-tp 
study the conductance spectrum for the spin current .t3 
Before presenting the results for the conductance, we 
give a physical analysis of the effect of exchange field on 
the AR. In Fig. g, we schematically draw the spin-split 
energy band in the ferromagnet and the AR process in 
the continuum model. As shown in the figure, the in- 
cident electrons and the Andreev reflected holes occupy 
different spin bands. Thus the AR is sensitive to the rela- 
tive spin-dependent density of states at the Fermi energy 
Ep. In the normal metal (/iq = 0), the energy band is 
spin degenerate, the AR is thus not suppressed. How- 
ever, if the exchange field is sufficiently strong that there 
are no electrons occupying the spin-down band in the 
ferromagnet, the incident spin-up electrons have no spin- 
down electrons to drag in order to form Cooper pairs. As 
a consequence, the AR is completely depressed. For the 
numerical calculation, we take Ny = 625. In Figs. and 
^, the subgap conductance spectrum G versus the scaled 
energy E is plotted for the FS junction and the FD^y 
junction with various values of ho but without the inter- 
facial scattering. As can be seen, for both the FS and 
FDxy junctions, the averaged conductance at a given en- 
ergy E is suppressed due to the blocking of AR. However, 
because the effective energy gap for the d-wave pairing 
symmetry is momentum-dependent while that of s-wave 
pairing symmetry is a constant in the momentum space, 
the different conductance behaviors between the FS and 
FD^y junctions are exhibited. In the FS junction, be- 
fore the subgap conductance is completely suppressed, a 
flat zero-bias maximum always shows up in the conduc- 
tance spectrum. When the exchange fleld is sufSciently 
strong, the conductance is zero within the energy gap, 
and sharply increases to a flnite value as the bias crosses 
the gap edge. This is because, outside the energy gap, 
the normal conduction process becomes important. In 
the FFlj-y junction, as the exchange fleld becomes strong, 
a zero-bias conductance maximum gives way to a zero- 
bias conductance dip. The strikingly similarity between 
the lowest curve in Fig. ^ and the experimental measure- 
ment performed on the La2/3Bai/3Mn03/DyBa2Cu307 
junctionsu demonstrates that the high degree of spin po- 
larization in the doped lanthanum manganite compounds 
and the d-wave pairing symmetry of the high- Tc super- 
conductors are essential to explain the observed conduc- 
tance behavior. In addition, as is shown, the conductance 
spectrum in the ferromagnet-superconductor junction is 



symmetric to the zero bias, which has been ohseraed 
in many experiments on spin polarized transport .Elij'a A 
general proof of this property is given in the Appendix W. 
In Figs. ^ and ^, we plot the conductance spectrum 
for a variety values of exchange field with the barrier 
strength flxed at Uq ~ Q.2D in the FS junction and 
Uq = 0.625D in the FD^y junction. In this case, the 
overall conductance spectrum is also reduced by the in- 
crease of ho. In the FS junction, a gap- like structure 
is exhibited in the conductance, and the peak at the 
gap edge is remarkably depressed by the exchange fleld, 
which is consistent with the recent experimental obser- 
vations on ferromagnet-s-wave superconductor junctions 
where the degree of spin polarization is small and a small 
barrier scattering potential may still exist B In the FD^y 
junction, due to the existence of midgap states at the 
interface, a sharp zero-bias conductance peak (ZBCP) 
shows up. The amplitude of this conductance peak is 
strongly suppressed by the exchange field. Meanwhile, 
as the exchange field becomes much stronger, the highly 
suppressed ZBCP is split. Physically, the suppression of 
the d-wave order parameter near the interface allows the 
ferromagnetic effect to penetrate into the superconduc- 
tor side through the tunneling of electrons, which leads 
to a small imbalance of the local occupation of electron 
with different spin direction so that a small magnetiza- 
tion at the d-wave superconductor side is induced. The 
small magnetization in tmn causes the shift of the en- 
ergy of the midgap states and the conductance peak is 
split. This splitting depends on the transparency of the 
interface. For a very strong barrier, the splitting is al- 
most unobservable. The splitting of the ZBCP by the 
exchange interaction can also be realized by the applica- 
tion of a magnetic field. If an in-plane magnetic field B is 
applied parallel to the interface, the orbital coupling be- 
tween electrons and the magnetic field can be neglected 
and only the Zeeman coupling ^^bB (fiB is the Bohr 
magneton) is present. Unlike the ferromagnetic effect on 
the electronic structures in the superconducting region 
near the interface, which is essentially of dynamic ori- 
gin, the Zeeman coupling is purely a local interaction. 
In this situation, the electron energy globely shifts to 
E ± ij^bB so that the energy of the midgap states shifts 
IJ-bB. Consequently, as shown in Fig. O, the ZBCP in 
the normal-metal-d2,2_ 3,2 -wave superconductor junction 
with {110}-oriented interface (From now on we call it 
the NDxy junction) can be readily split. The range of 
splitting is just 2^bB. 



B. The effects of spin- flip interfacial scattering 

In the preceding treatment, the spin-flip interfacial 
scattering effects are ignored. In case of the junctions 
with the ferromagnet involved, this type of scattering 
may be important. To study the effect, we introduce a 
new term into the Hamiltonian 



'Hsp 



E 



ut 



(4.11) 



where Uf^^, = UidonSaa' is assumed to be nonzero at 
the interface (layer index n = 0). In the spin-space, the 
spin-flip scattering term is represented by 



U'P = 



Ui 
Ui 



(4.12) 



The spin-nonflip term has been represented by Ui in the 
total Hamiltonian, which is the diagonal elements in the 
spin space. In the normal state junction with the spin- 
flip scattering, a beam of incident electrons with spin 
a will be reflected as electrons with the same spin and 
opposite spin. In the junction made up of the supercon- 
ductors, one can expect that, when a beam of electrons 
with spin a incident from the normal metal or ferromag- 
net, the spin-flip scattering will also lead to the normally 
reflected electrons with the opposite spin g_ and Andreev 
reflected holes with the same spin cr, in addition to the 
reflected electrons with the same spin a and holes with 
the opposite spin g_. By working with the coupled 4x4 
BdG matrix equations, we can generalize the previous 
transfer matrix technique to obtain the normal reflection 
amplitudes, r^^aa, fe^aa-, and the Andreev reflection am- 
plitudes, r^.aa, rh,crcr- Correspondingly, the conductance 
is generalized to be 



G 



hN^ 



/ ^ [J- 1 ^h^aa I ^h^aa ^e,crcr ^e^acrl ; 



shows up. (see Fig. |T^)(c)) As the spin-flip scattering 
strength is comparable to the spin-nonflip part, the com- 
plementary behavior in the variation between Re,aa and 
Re.aa and weak bias-dependence of Rh.aa (highly sup- 
pressed within the gap) and Rh,aa leads to an almost 
constant conductance spectrum, (see Fig. |2|(d)) If the 
spin-flip scattering is much stronger than the spin-nonflip 
part (Ui = ID), the induction of a peak at finite bias in 
both Rh.aa and Rh,aa causes a finite-bias conductance 
peak, (see Fig. |l^(e)) Whether such a extremely strong 
spin-flip scattering compared with the spin-nonflip scat- 
tering exits experimentally is unclear and we will not 
discuss this extreme limit further. As shown in Fig. M, 
the ZBCP can be completely depressed in the averaged 
conductance spectrum of a NDj;y junction [Uo = 0.625Z3) 
by the spin-flip scattering. Figure nj plots the averaged 
conductance for various values of Ui in a VY>xy junction 
with C/o = 0.625Z? and ho = 0.4751?. Since the spin- 
flip interfacial scattering tends to spoil the pre-oriented 
spin direction of conduction electrons incident from the 
ferromagnet, it seriously influences the conductance spec- 
trum. In particular, the splitting of the ZBCP induced 
by the exchange fleld is washed out in the presence of 
a strong spin-flip interfacial scattering so that the con- 
ductance spectrum becomes completely structureless. In 
addition, as shown in Figure |lj, the suppression of the 
conductance at the region away from zero bias by the 
exchange field is reduced by a strong spin-flip interfacial 
scattering. 



where Rh^ap = [sin((7,i^o.a)/sin((je^c,a)]|rft,o 

[sin{qh.,ja)/ Sa\{qe.aa)\\rh,aa?, Re.aa = 



(4.13) 



and 



Re. 



= [sin((je,ffa)/sin(ge,o-a)]ke 



with g's the wave 



vectors associated with different types of electrons and 
holes. In Fig. |l^, we plot the conductance spectrum and 
the reflection coefficients for various values of the spin- 
flip interfacial scattering strength Ui in a NDajj, junc- 
tion with the spin-nonflip interfacial scattering strength 
Uq = 0.625Z3. The transverse momentum is taken to be 
ky = 7r/3V2a. In the absence of the spin-flip scattering, 
i.e., Ui — 0, the coefficients Re,aa and Rh,aa vanish, and 
Rh,aa decreases while Re^a increases monotonically with 
the bias within the effective energy gap |Ak| = 2Arfo- 
Especially, Rh,aa = 1 and Re,atT = at i? = 0, which ac- 
counts for the appearance of the ZBCP. (see Fig. p^(a)) 
The effects of the spin-flip scattering on the conductance 
spectrum depends on its strength in detail. For a rela- 
tive small value of the spin-flip scattering (C/i — 0.125D), 
Re.aa and Rh,aa dccrcascs monotonically with the bias. 
Re,aa IS flnitc and Rh,aa is depressed aX E = Q, which 
leads to a suppressed ZBCP. (see Fig. |2|(b)) As the spin- 
flip scattering strength is further increased (C/i = 0.251?), 
Rh,aa varies nonmonotonically, it flrst increases to reach 
a maximum and then decreases with the bias. In addi- 
tion, the zero-bias Re,aa and Re,aa are enhanced. Con- 
sequently, a fiat conductance maximum at a finite bias 



V. LOCAL MAGNETIC MOMENT 

It has been predictecO that besides the ZBCP in the 
quasiparticle tunneling, one of the other consequences of 
midgap states is the possibility of a sizable magnetic mo- 
ment at the {110} surface of the d^2_^2--wai,Ye supercon- 
ductor. Actually, the splitting of the ZBCP in the FD^y 
junction with a strong exchange field or in the NDa;j, junc- 
tion with an in-plane magnetic field haS|-supported this 
prediction. In a recent theoretical work,E3 a formal ex- 
pression for the magnetic moment has been given, but 
a serious calculation of this quantity has not been done. 
In this section, we give a detailed analysis of the local 
magnetic moment (LMM). 

The average electron density for each spin direction is 
given by 

(nia) = Y.{\ul\'fiE,) + KHl - fiE,)]} , (5.1) 

V 

where ]{E) = [l + exp(_E/r)]^^ is the Fermi distribution 
function. The local magnetic moment can be defined as 



m = -~^iB{ni^ -"-ii) 



(5.2) 



For the FD^y junction at temperature T = O.OSA^o, the 
LMM at the distance x = a/^/2 away from the inter- 
face in the superconducting region is found to be: When 



ho = 0.5D, m{x = a/V^) = O.OlAfiB for Uq = 0.6251), 
and 0.008/z_B for t^o = 2.5D. Corresponding to the split- 
ting of the ZBCP in the FD^^j, junction, the exchange- 
field induced local magnetization at the surface of d- 
wave superconductor is sensitive to the interfacial bar- 
rier strength. In Fig. |l^, we plot the magnetic-field de- 
pendence of the LMM at a; = a/y/2 for different tem- 
peratures. At low fields, the LMM varies linearly with 
the field. The slope increases with the decreased tem- 
perature. At higher fields (about six times of the tem- 
perature) so that the width of the midgap peak in den- 
sity of states is surpassed, the LMM begins to saturate. 
In contrast, for the FS junction or the normal-metal-s- 
wave superconductor junction with a Zceman coupling, 
we find that the LMM is almost zero. In Fig. |l^, we 
plot the spatial variation of the LMM into the super- 
conducting region of the NDj,j, junction with Uq — 2.5D 
and iibB = 0.4Ado. The temperature T = O.OSArfo- It 
is shown that the LMM has nonzero values only at the 
sublattice x — [2n -I- l)a/y/2 with n being non-negative 
integer, and it decays into the bulk of the superconductor. 
This interesting behavior directly reflects the existence of 
zero-energy peak at x = (2n -I- l)a/-\/2 and the absence 
at a; = 2na/-\/2 as the chemical potential /i = 0H3 This 
feature is special to the lattice model. 



VI. CONCLUSIONS 

In conclusion, we have presented a unified theory for 
the proximity effect, quasiparticle transport, and lo- 
cal magnetic moment in the FDj;^ junctions by solving 
the Bogoliubov-de Gcnnes equations within an extended 
Hubbard model. As a comparison, the calculations are 
also made for the FS junctions. The energy gap ap- 
pearing in the BdG equations have been determined self- 
consistently by using exact diagonalization technique. It 
is found that the proximity induced order parameter os- 
cillates in the ferromagnetic region but is almost un- 
changed in the superconducting region by the exchange 
field. The modulation period of the proximity induced 
order parameter is shortened by the exchange field but 
the oscillation amplitude is decreased with the interfa- 
cial scattering. Once the superconducting energy gap 
for various interfacial scattering potentials is determined 
self-consistently, a transfer matrix method has been pro- 
posed to calculate the the subgap conductance within a 
scattering approach. We find that the subgap conduc- 
tance is suppressed by the spin splitting of the energy 
band in the ferromagnet. For a ballistic FDa;^ junction, a 
conductance dip is exhibited with strong exchange fields. 
In the presence of interfacial scattering, the ZBCP is 
split by the strong exchange field. The degree of this 
splitting depends on the barrier strength. In contrast, 
the ZBCP can be split very easily by an in-plane mag- 
netic field due to the local nature of the Zeeman coupling 
and the range of splitting is independent of the barrier 



strength. In addition, we also show that the spin-flip in- 
terfacial scattering can seriously influence the quasipar- 
ticle transport properties. As one of the consequences of 
the midgap states, a sizable local magnetic moment in 
the FDa;j, junction or in the NDa;^ junction in the pres- 
ence of the Zeeman coupling has been found. Inspired by 
the observation of the zero-bias conductance dip in the 
ballistic ferromagnet-high- Tc superconductor junctions, 
we believe that the other interesting behaviors predicted 
in this paper are also experimentally accessible. 
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APPENDIX A: SYMMETRY PROPERTY OF 
THE CONDUCTANCE SPECTRUM 

T he e igenstates of the BdG equations give n by 
Eq. (2.6) exist in pairs: Each eigenst ate fr om Eq. ( 2.6a| ) 
is related to a counterpart from Eq. ( 2.6b ) as: 






-E 



(Al) 



This mirror image property makes the conductance spec- 
trum is symmetric about the zero of bias. Suppose 
we have a incident, and outgoing waves as solutions to 
Eq. (2.6a) with energy E 



and 



p-iqhix 






The amplitudes of the incoming and outgoing waves are 
then connected by the scattering matrix: 






^out 
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From the mirror image property of the eigcnfunctions, 
the incident and outgoing waves at energy — _E must be 
given by 



and 






-h* p-^lhi'^ 



These wavefunctions is related by the scattering matrix 
at energy —E: 



'''h],ei rh],h] 
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The above equation can be rewritten as: 



(A3) 



fh],h] 






(A4) 



Comparing Eq. A4 with Eq. A2, we find the relation 



M,eT(i?) = KTJ.T(-^)]% 

rhu^{E) = [-r,^ji^{-E)r 
rHiM{E)^[rei.ei{-E)]* . 



(A5) 



The symmetry properties of the reflection amplitudes 
yields the identity: 

G„{E)^G„{-E). (A6) 

Therefore, we have extended the symmetry propertyE3 to 
the magnetic case. 
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FIG. 10. The differential conductance spectrum for var- 
ious values of exchange field in the FDxy junction with 
Uo = 0.625D. 



FIG. 1. Schematic geometry of the ferromagnet-s-wave 
superconductor junction with a {100}-oriented interface (a) 
and the ferromagnet-d-wave superconductor junction with a 
{110}-oriented interface defined on a two-dimensional lattice. 
The filled circles represent the ionic positions in the ferromag- 
net and the filled squares represent the ionic positions in the 
superconductor. The interface layer is represented by empty 
circles. Xa.6 are crystalline axes. The propfiles of s-wave and 
d-wave order parameters are also shown, respectively. 



FIG. 2. The spatial variation of order parameter for various 
values of exchange field in the FS junction without the inter- 
facial scattering potential. Here D = 8t is the band width. 



FIG. 3. The spatial variation of order parameter for various 
values of exchange field in the FDxy junction without the 
interfacial scattering potential. 



FIG. 4. The spatial variation of order parameter for various 
values of interfacial scattering potential in the FS junction 
with the exchange field ho = 0.125D. 



FIG. 5. The spatial variation of order parameter for various 
values of interfacial scattering potential in the FDxy junction 
with ho = 0.125D. 



FIG. 11. The differential conductance spectrum for vari- 
ous values of Zeeman coupling /is in the ND^jy junction with 
Uo = 0.625D. 



FIG. 12. The transverse-momentum-dependent differen- 
tial conductance spectrum and the reflection coefficients in 
a NDa,y junction for various values of spin-fiip scattering 
strength Ui = (a), 0.125D (b), 0.25D (c), 0.625D (d), and 
ID (e). The spin-nonfiip scattering strength Uo = 0.625-D. 
The transverse momentum ky — ■K/?)\/2a. The effective en- 
ergy gap for this momentum is |Ak| = 2Ado- 



FIG. 13. The averaged differential conductance spectrum 
in a NDajy junction for various values of spin-fiip scattering 
strength. Uo = 0.62573. 



FIG. 14. The averaged differential conductance spectrum 
in a FDxy junction for various values of spin-flip scattering 
strength. Uo = 0.625LI and ho = 0A75D. 



FIG. 15. The magnetic-field dependence of the local mag- 
netic moment at the distance a/\/2 away from the interface 
in superconducting region of the ND^,;, junction at different 
temperatures. The interfacial scattering potential Uo ~ 2.5D. 



FIG. 6. A schematic drawing of (a) the spin-split energy 
band in the ferromagnet within the continuum model and (b) 
the Andreev refiection process at the interface between the 
ferromagnet and superconductor: A beam of spin-up elec- 
trons incident with angle 6m and energy within the gap are 
normally reflected as spin-up electrons and Andreev reflected 
as spin-down holes. The Andreev reflection angle 6a is re- 
lated to 9n by the conservation of the momentum component 
parallel to the interface. The thick solid line represents the 
interfacial scattering layer. Also shown are the d-wave order 
parameter proflle and the angle a of the crystalline orienta- 
tion with respect to the interface. Ef is the Fermi energy. 



FIG. 16. The spatial variation of the local magnetic mo- 
ment into the superconducting region of the ND^jj; junction 
at T = O.OSAdo and hbB — 0.4Ado- The interfacial scatter- 
ing potential Uo = 2.5D. 



FIG. 7. The differential conductance spectrum for various 
values of exchange field in the FS junction without the inter- 
facial scattering potential. 



FIG. 8. The differential conductance spectrum for various 
values of exchange field in the FDxy junction without the 
interfacial scattering potential. 



FIG. 9. The differential conductance spectrum for various 
values of exchange field in the FS junction with Uo ~ 0.2D. 
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